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Abstract 

We study the scattering of an electron from a one dimensional inverted Gaussian atomic 
potential in the presence of strong time periodic electric fields. Using Floquet theory, we 
construct the Floquet Scattering matrix in the Kramers-Henneberger frame. We compute 
the transmission coefficients as a function of electron incident energy and find that they 
display asymmetric Fano resonances due to the electron interaction with the driving field. We 
find that the Fano resonances are associated with zero-pole pairs of the Floquet Scattering 
matrix in the complex energy plane. Another way we "probe" the complex spectrum of the 
system is by computing the Wigner-Smith delay times. Finally we find that the eigenphases 
of the Floquet Scattering matrix undergo a number of "avoided crossings" as a function of 
electron Floquet energy, and this number increases with increasing strength of the driving 
field. These "avoided crossings" appear to be quantum manifestations of the destruction of 
the constants of motion and the onset of chaos in classical phase space. 



PACS numbers: 32.80.-t, 34.50.-s, 03.65.Nk, 05.45.Mt 



I. INTRODUCTION 

The development of ultra-high intensity lasers has led to the study of atoms in external 
time periodic electric fields that are comparable in strength to the electric fields produced by 
the atomic nucleus. One of the most interesting phenomena observed in these time periodic 
systems is the stabilization with increasing laser intensity that was predicted theoretically 
|T|-^ and has been verified experimentally [§,0. 

In previous studies, one dimensional atomic potentials have been used to predict sev- 
eral phenomena in the theory of laser-atom interactions at high laser intensities. Many of 
these studies were carried out in the context of Floquet theory formulated in the Kramers- 
Henneberger (K-H) frame of reference , which oscillates with a free electron in the time 
periodic field. Gavrila and Kaminski developed a nonperturbative method to study elec- 
tron scattering in the presence of strong time-periodic electric fields. Using three dimensional 
models, Dimou and Faisal as well as Collins and Csanak have studied resonances in 
laser assisted scattering. Bhatt, Piraux and Burnett |P in their work on electron scattering 
from a polarization potential in the presence of strong monochromatic light, argued the 
appearance of new light-induced quasibound states (resonances) as the field strength is in- 
creased. The same phenomenon was later also observed by Bardsley and Comella |TI]] and 



Yao and Chu [|TT| who used the complex coordinate scahng transformation to compute the 
complex quasibound states in their study of photodetachment from a one-dimensional Gaus- 
sian potential. In addition, the same atomic potential was used by Marinescu and Gavrila 
|l^ to compare the predictions of the full Floquet theory with those of the high-frequency 



Floquet theory (HFFT), using resonance (Siegert) boundary conditions. (The HFFT theory 
is a version of the Floquet approach adapted to treat the high frequency limit.) Recently, 
Timberlake and Reichl using the inverted Gaussian potential, studied the phase-space 
picture of resonance creation and they showed that the light-induced quasibound states are 
scarred on unstable periodic orbits of the classical motion. 

In this paper, we study the scattering of an electron from a short range atomic potential 
in the presence of a strong time periodic electric field. The atomic potential we consider is 
the one-dimensional inverted Gaussian potential, a model that has already offered interesting 
insights into different aspects of the laser-atom interactions [p!0|-[l2[|. Our goal is to construct 
the Floquet Scattering matrix (S-matrix) using the full Floquet theory formulated in the 
K-H frame for a strongly driven atomic system. A Floquet Scattering matrix has been 
constructed by Li and Reichl [|1^] for periodically driven mesoscopic systems. The Floquet 
S'-matrix connects the outgoing propagating modes to the incoming propagating modes and 
is a unitary matrix which conserves probability. We construct the Floquet S-matrix in the 
K-H frame where we can define asymptotic states. 

In section II, we construct the Floquet S-matrix in the K-H frame. In section HI. A we 
compute the transmission coefficients and the poles of the Floquet S'-matrix and find the 
quasibound states of the atomic system. In section III.B we compute the Wigner-Smith delay 
times of the scattered electron as a function of the electron incident energy and show that the 
Wigner-Smith delay times of the scattered electron due to the presence of the quasibound 
states are of the same order of magnitude as the lifetimes obtained from the poles of the 
Floquet S- matrix. We finish, in section III.C, with a quite interesting observation. When 
plotting the eigenphases of the unitary Floquet S-matrix as a function of the electron's 
Floquet energy we find that at certain energies the eigenphases undergo "avoided crossings" 
that change the eigenphases character completely. We find that the number of "avoided 
crossings" increases with increasing strength of the time periodic electric field. The "avoided 
crossings" observed as the strength of the driving field is increased appear to be quantum 
manifestations of the destruction of the KAM surfaces and the onset of chaos in the classical 
phase space. 



II. THE FLOQUET S-MATRIX 
A. The model 

We study the scattering of an electron in the presence of a strong electric field and a 
short-range atomic potential. The electric field E{t) = Eq sin(u;t) (T = — is the period of 
the field) is treated within the dipole approximation as a monochromatic infinite plane wave 
linearly polarized along the direction of the incident electron. The Schrodinger equation, in 
one space dimension x, that describes the dynamics of the system is in atomic units (a.u.) 
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= (^(-^^ - iMt)r + v{x)^ t), (1) 

where V(x) is the inverted Gaussian potentiah 

V{x) = -Voe-^^/'^\ (2) 



and q is the particle charge which for the electron is q = —1 a.u.. The electric field is, 

at ' 



E{t) = — ^^f^ , where A(t) is the vector potential and is given by 



E 

A{t) = —cosiut). (3) 

CO 

We use atomic units {e = h = m = 1) throughout this paper, except when otherwise 
indicated. 

To construct the Floquet S'- matrix of the system, we transform to the K-H frame 
In the K-H frame there are well defined asymptotic regions and the boundary conditions are 
expressed in terms of free electron waves. To obtain the wave function in the K-H frame we 
introduce the unitary transformation ^,0] 

^x,t) = UiU2^{x,t), (4) 

where 

U, = e*/-^^^(*')'^*' and U, = ^-^^^''>^'' ■ (5) 



Ui is a phase transformation to remove the term from Eq.(|lD while U2 is a space trans- 
lation transformation to the K-H frame. In the K-H frame, the wave function satisfies the 
following Schrodinger equation 



d<^{x,t) ( 1 d- 



2 



dt \ 2 



+ V{x + a{t))]^{x,t), (6) 



where a{t) is the classical displacement of a free electron from its center of oscillation in the 
time periodic field E{t), and is given by 

a{t) = -g/!^ A{t')dt' = ao sin(tut) with ao = ^ . (7) 

In Eq.(^, the potential is a periodic function of time, that is V{x+a{t)) = V{x+a{t+T)). 
Thus, according to the Floquet theorem |T^, Eq.(^) has solutions of the form 

Mx,t)=e-'''Mx,t), (8) 

where £ is the Floquet energy, £ G [0,ti;), and (j)s{x,t) is a periodic function of time, 
(f)£{x, t) = t + T). Taking the Fourier expansion of (psi^x, t) we obtain 

+00 

^£ix,t)=e-''' Y: ct>n{x)e-''^\ (9) 
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where n indicates the Floquet channel. Note that (f)n{x) is also S dependent but we omit the 
S subscript to simplify notation. The energy E of an incident electron in the K-H frame is 
related to the Floquet energy through the expression E = S + nuj. Next, we Fourier analyze 
the potential 

+00 

V{x + a{t))= K(ao;a;)e-*""'*, (10) 

n=— 00 

where the Fourier components for the inverted Gaussian potential, Eq.(0), can be written 

as 

Vn{ao;x) = — r V{x + a{t))e'''^'d{ujt) = ~Vq- T cos(ncjt)e-(^+"»™^('"*))'/'''c/(cut), (11) 
271 Jo vr JO 

see Fig.(|I]). To be able to construct the Floquet S-matrix, the Fourier components \4(ao; 
must be smooth functions in the one space dimension x in the K-H frame. This is indeed the 
case for the inverted Gaussian potential. Eg . (|Tl|) . Note that in the K-H frame the potential 
oscillates back and forth along the x-axis (laterally) with the period of the external field. 

From Eq.(|lT]), we see that the components V^(ao; x) of the atomic potential in the limit 
X — > ±00 tend to zero faster than 1/x and we can thus divide the one space dimension 
X in three regions: the asymptotic regions I, x G [xo,oo), and HI, x G (— 00,— xq], where 
the potential is asymptotically zero; and the scattering region II, X G [ — Xq , Xq 1, where the 
potential V{x + a(t)) is not zero, see Fig.(||). In the rest of this paper, for brevity, we refer 
to the potential in asymptotic regions I and HI as being zero instead of asymptotically zero. 
The choice of xq depends on the value of the parameter ao. The larger ao is, the further out 
we have to define the asymptotic regions I and HI. 



B. Floquet solution in the scattering region II 

Substituting Eqs.(|D and ([l0| ) into Eq.(|D we obtain an infinite system of coupled differ- 
ential equations [Q for the Floquet components (pn{x) 

1 +°° 
- 2~^^^^^) + [^o(ao; x) -{£ + nuj)\(t)n{x) + ^ K-z(ao; 2;)0z(x) = 0. (12) 

/ — — 00 

Next, we truncate to a finite number of Floquet channels and take Ue and Hp to be the lower 
and upper limit of the Floquet channels considered. That is, n = —Ue, ...,0, ...,np and the 
total number of Floquet channels is given hj N = Ue + Up + 1. After truncating, Eq.(|12|) 
can be cast in the following matrix form, 

I^2^'\^) = M{x)4>'\x), (13) 

where / is the unit N x N matrix, (f)^\x) is the x 1 matrix with elements 0^^(x) = (f)n{x) 
and M{x) is an iV x iV matrix with elements 

Mn,iix) = 2(K-K«o; x) - 5n,ii£ + nu)), (14) 
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where Sn,i is the Kronecker delta and n, / = —He, •••,0, ...,np. The general solution of the sec- 
ond order N coupled differential equations, Eq. (p!3D , can be written as a linear combination 
of 2N linearly independent columns with j = 1, ...,2N, as follows 

= X^^\x)C + X^^\x)D (15) 

where X^^\x) and X^'^\x) are N x N matrices whose elements are functions of the one 
space dimension x and C, D are constant N xl matrices. Each of the linearly independent 
columns Xji^) satisfies Eq.(|TB|) 

^-^^A^) = M{x)Xj{.x) -^Xnjix) = MnAx)Xl,j{x), (16) 

l=—ne 

where n = —He, 0, np and j = 1, 2N. The functions Xn,j{x) can be found analytically 
if the matrix elements of M{x) are constant. From Eq.(|l^), it follows that every channel 
function 0^^(x) can be written as a linear combination of 2N functions Xn,j{x) and thus, the 
wavefunction in the scattering region II is given by 

Up N 
n=—ns 771=1 



C. Floquet solution in the asymptotic regions 

In the asymptotic regions I and III the potential V{x + a{t)) is zero. Thus, we can 
consider as our boundary conditions a superposition of incoming and outgoing free electron 
waves in the truncated Floquet channels that are incident from both sides of the scattering 
region: 

"p '^P piknX "p p—iknX 

\ ^ iIf„,\^—iCt^—imjjt \ ^ uout „—i£t „—inu)t i \ ^ uin- „—iCt „—inu>t 



n=-rie n=-ne V"^n n=—ne V "^n 



-ik„x piknX 



(18) 



(19) 

where 6^", a^" and 6°"*, a°"* are the probability amphtudes of the incoming and outgoing 
electron waves, respectively, that are incident in the nth Floquet channel with energy E = 
S+nuj, see Fig.(^. Propagating modes are incident on the Floquet channels n = 0, rip and 

have wavevectors fc„ = ^J2{S + nu), while evanescent modes occupy the Floquet channels 

n = — ne,...,— 1 and have imaginary wavevectors kn = i\j2\£ + nijj\. The current density 
of the evanescent modes is zero. We note that the terms propagating/evanescent modes 
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correspond to what some authors refer to as open/closed channels, respectively. For the 
Floquet S-matrix to be unitary we need to normalize the current density of the propagating 
modes. To do so, we have introduced the constants in the wavefunction in Eqs.(|18[) 

and ([l^) . To simplify notation in Eqs. (|T^) and (|T^) , we introduce the constants for 
the evanescent modes as well, even though they have zero current density. 

It is important to note, once again, that the reason we choose to work in the K-H frame is 
that in this frame we can define asymptotic regions where the potential is zero and thus, the 
Floquet channels are not coupled, in contrast with the scattering region, as we have already 
shown. The existence of the asymptotic regions guarantees that probability is conserved 
in the truncated number of Floquet channels and thus the Floquet S-matrix is a unitary 
matrix. 



D. Floquet S-matrix 



The Floquet S-matrix connects the outgoing propagating modes with the incoming prop- 
agating modes, and in this section we show how to construct it. As we show in what follows, 
the Floquet S^-matrix connects channels with energies that differ by an integer multiple of 
uj, while in the usual time-independent scattering theory the S-matrix connects channels 
with the same energy. The reason is that the Floquet S-matrix describes a time- dependent 
process and thus the energy of the incident electron is not conserved. However, because the 
Hamiltonian is time periodic, according to Floquet theory [|T5|, the Floquet energy £ defined 
modulo u; is a conserved quantity. 

The wavefunction and its first spatial derivative must be continuous at the boundaries 
of the asymptotic regions ±xo. At x = xq these conditions lead to 

^r*^^ + n- = {Xn,mMCm + Xn,N+m{XQ)dm) , (20) 



iknb. 



out 



-iknXo 



N 

E 

m=l 



\'\.n,r, 



{Xo)Cm + Xn,N-\ 



{Xo)dm), 



(21) 



where x'n,J^) = ^^^^ and Xn,N+mi^) = """^'^ , ^^il^ at a; = -Xo they lead to 



,out 



xo 



xo 



+ al 



N 

E 

m=l 



(Xn,m( ■^o)(^m ~l~ Xn,A''+m( '^o)'^m)) 



(22) 



- ikna°^^^= + ifcn<" rrr = E (Xn,m(-2;o)Cm + Xn,N+m{-Xo)dm) ■ (23) 
V l^n V l^n m=l 

Due to the connection conditions Eqs.(|g), (|TD, (H) and(|2|) only 2N out of the QN 
coefficients are arbitrary and we choose those to be the incoming probability amplitudes 
a™ and 6™. In Eqs.(^, (^, (|2^) and(|^), the probability amplitudes a*" and of the 
evanescent modes are zero because of the unbounded character of the exponentials they 
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multiply in the asymptotic regions I and III. That is, 6^" = a^" = 0, for n = —rie, — 1. We 
now introduce the N x N matrices 

iK±)n,l = e^'^-''°5n,U {K'±)n,l = ±lkne^"'-'^' Sn,l, U, I = -He, 0, .., Up, (24) 



T _ j 0, if 7^ / and if = / = —rig, —1, 



1, if n = Z = 0, Up 



J^n,i= 7h^n,h n,/ = -ne,...,0,...,np. (26) 



and 



Also ^ X«(±xo), X'i^ ^ X(2)(±xo), X«l = X«'(±xo), X^^^^ = X(2)'(±xo), 

where X^^^ (x) and X*-^-* (x) are the derivatives of X^^\x) and X^'^\x) with respect to the 
one space dimension x. We also introduce the x 1 matrices A°"* = a°"*, = 6°"*, 

A^" = aj,", S^" = fej;^. Next, we write Eqs.(^), (^), and (^ in matrix form as follows 

A^K+S""* + ATK^ JS^" = X^+^ C + X^D (27) 
J\fK\B°''' + J\fK'_JB'' = X(i)^C + (28) 
A^K+A""* + J^K_JA''' = X'^I^C + xL^^i:) (29) 

- AfK'^A""' - AfK'_JA"' = X(^)'_C + X^^^Id (30) 

After algebra given in Appendix A we find the Floquet S-matrix, that connects the outgoing 
probability amplitudes of the propagating modes to the incoming probability amplitudes of 
the propagating modes, to be 

(•^p-p "^pp^pp ^ pp ^pp^ pp 
pp ''pp-'^ pp pp ' pp-'^ pp 






(31) 

where the Up + 1 x Up + 1 {rip + 1 is the number of the propagating modes) matrices r^^, 
Vpp, t'pp and tpp defined in Eg s . (1461) and (^7|) of Appendix A connect propagating modes 



to propagating modes and contain the evanescent mode effect as is shown in Appendix A. 
Also the rip + 1 X 1 matrices A^" A""**, B^^ and have elements the amplitudes of the 
propagating modes and are defined in Eq. (^) of Appendix A, and the + 1 x np + 1 matrix 
Afpp has elements the normalization constants of the propagating modes and is defined in 
Eq.(^8D of Appendix A. The matrices R' , R, T' and T have dimensions + 1 x np + 1 
and their elements are given in terms of the elements of the r^^, Vpp, t'^^ and tpp matrices as 
follows 
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^^n\n V kn \' pp)^^"^^ JMi',n — \J kn \' PP)n',n, n' ,n ~ \J \ PP^^'"-^ -Ln',n — \ k„ \''PP)n',n ! 

(32) 

with n', n = 0, Hp. 

In Appendix B we show how to obtain numerically the matrices Vpp and tpp. The Floquet 
S-matrix has dimensions 2(np+l) x2(np+l), see Eg. (|HTD , and is determined by the reflection 
and transmission amplitudes, R'^, T'^i „ and T„/ „, of the propagating modes. The 

elements |-Rn',nP and |T„/^„p are the reflection and transmission coefficients respectively for 
an electron wave incident on the propagating channel n from the right that gets scattered 
to the propagating channel n', while the elements and |T^/„P are the reflection 

and transmission coefficients respectively for an electron wave incident on the propagating 
channel n from the left that gets scattered to the propagating channel n' . 

In this section we have shown how to construct the Floquet S'-matrix in the K-H frame. 
The reason we work in the K-H frame is that we can define asymptotic regions where the 
wavefunction is a superposition of free electron waves. That guarantees that the truncated 
Floquet S-matrix is a unitary matrix, that is, the following condition is satisfied 

Up 

[|-Rn',nP + l^n'.nP] = 1, (33) 

ra'=0 

for every incident propagating mode n = 0,...,np. The above condition is a statement of 
conservation of probability. Also, the Floquet S-matrix we construct in the K-H frame is 
isospectral with the corresponding matrix in the Lab frame since a unitary transformation is 
used to transform from the Lab to the K-H frame, see section II. A. Finally, the criterion we 
use to successfully truncate to N Floquet channels is that an electron wave incident on the 
last propagating Floquet channel n = Up is not affected by the scattering potential. That is, 
the transmission coefficient {Tn^^npl"^ should be equal to one as a function of electron incident 
energy E { E = S + Upuo) as we discuss in more detail in section III.B. 



E. Symmetries of the Floquet S-matrix 

The Hamiltonian of the scattering model we consider in the K-H frame, Eq.(|^), is in- 
variant under the transformation x —x and t ^ t + T/2, which is known as Generalized 
Parity. Thus, H{x,t) = H{—x,t + T/2) and therefore $(— x,t + T/2) is also a solution of 
Eq.(^. Applying the transformation x —x and t ^ t + T/2 to Eqs.(|T^) and (|T^) it is 
easy to show that the Floquet S-matrix has the following symmetry: 

K',n = i?„',n(-l)"'-", K>,n = T„',„(-l)"'-", with n',n = 0, ...,np . (34) 

Thus, if we know the reflection/transmission amplitudes, Rn',n/Tn',n for electron waves in- 
cident from the right using Eqs. (|3^) we can find the reflection/transmission amplitudes, 
i?^/ „/T^/ „ for electron waves incident from the left and vise versa. 
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III. RESULTS 



In this section, the calculations are performed with the values Vq = 0.27035 a.u. and 6 = 
2 a.u. assigned to the parameters of the inverted Gaussian potential. For these parameters 
the Gaussian potential supports only one bound state of energy Ef, = —0.1327 a.u. in the 
field-free case. The parameters Vq and 6 were chosen so as to describe the behavior of a 
one-dimensional model negative chlorine ion, Cl~ , in the presence of a laser field, and are the 
same as considered in 



11,12,16 



The frequency of the time periodic field is taken constant 
and equal to u = 0.236 a.u. for all our calculations. For these values of the parameters Vq, 
6 and u the inverted Gaussian potential has been shown to exhibit stabilization in|,[T2l . 



A. Transmission resonances 

In this section, we compute the transmission coefficient T„ and the total transmission 
coefficient Ttot,n as a function of the electron incident energy E, where 

r„(E) = \Tn,n?, TtotAE) = EnLo \Tn' ,n? , E G [uU , {u + 1)^) (35) 

with n = 0, 1. Thus, we consider an electron wave incident from the right with energy 
E G [0, 2uj) and compute the transmission coefficients. Keeping the frequency of the time 
periodic field constant, u = 0.236 a.u., and varying the strength of the driving field, oq, we 
plot the transmission coefficients r„ and T^j „ in Figs.(0), and (|^) for uq equal to 0.5, 
2.25 and 5.25, respectively. The frequency of the driving field u = 0.236 a.u. is chosen so 
that it is larger than the binding energy of the ground state \Ef,\ = \ — 0.1327| a.u. in the 
field-free case. The Floquet channels we retain to obtain the numerical results presented 
in sections III. A and III.B are n = — 6, 0, 6 for ao = 0.5, n = — 12, 0, 12 for 
ao = 2.25 and n = —19, 0, 19 for ao = 5.25, for reasons we discuss in detail at the end 
of section III.B. 

The transmission coefficients r„ and Ttot,n display sharp asymmetric resonances, as a 
function of electron incident energy E, that involve a dip or a transmission peak/dip as 
is shown in Figs.(0), (^) and @. These asymmetric resonances are due to the interaction 
of the incident electron wave with the laterally oscillating potential in the K-H frame and 
are the so-called Fano |T7|JT8[] type transmission resonances that are known to occur when a 



bound state is coupled to a continuum of states. This is indeed the case for the scattering 
model we consider, where the bound state of the inverted Gaussian potential is coupled to a 
continuum of states through the time periodic electric field. Note in Figs.@, (§) and (|) that 
the difference between the transmission coefficient r„ and Ttot,n becomes more prominent 
with increasing oq. The reason is that as ao is increased more Floquet channels interfere with 
the incident electron wave and significantly contribute to the total transmission coefficient. 
A comparison of Fig. (|) with Figs. (Q) and reveals that as the driving field is increased 
the higher order resonances, for E > uj, become stronger. 

We now focus on the transmission coefficient T„ and discuss how it "probes" the qua- 
sibound states of the system. For ao = 0.5, see Fig. (3a) the system has only one Fano 
transmission resonance which for small amplitude of the driving field is associated with 
the n = —1 localized Floquet evanescent mode which has its origin in the bound state of 
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the undriven system. When the strength of the driving field is increased, T„ has a second 
Fano transmission resonance at a higher incident energy, see Figs. (4a) and (5a). This sec- 
ond resonance appears for > 1, as was shown in |Tl|JT^, and it is thus a field induced 
resonance. 

The Fano-resonances, which are indicated by a dip or a transmission peak/dip in the 
coefficient Tn, correspond to quasibound states of the system that show up as poles of the 
Floquet S-matrix in the complex energy plane. In what follows, we compute the poles of 
Tn in the complex energy plane. Other elements of the Floquet S'-matrix have poles as 
well. As was noted in |]TU[ the asymmetric Fano line shape in T„ is associated with zero-pole 
pairs when plotting r„ in the complex energy plane. By a zero-pole pair we mean that 
every transmission zero of Tn along the real energy axis is associated with a pole of Tn on 



the lower half complex energy plane due to the unitarity of the Floquet S'-matrix For 
ao = 0.5 there is only one zero-pole pair associated with the single transmission resonance 
seen in r„, while for = 2.25 there is a zero-pole pair for each of the two resonances, see 
Figs. (4a) and (P). For small strengths of the driving field, = 0.5, the location of the pole 
on the lower half complex energy plane and of the zero on the real energy axis is the same, 
while there is a small difference for stronger fields, ao = 2.25. That is why, we can only 
approximately determine the real part of the quasibound states from the transmission zeros. 
From the poles in the complex energy plane, we find the real part of the quasibound states 
to be Re{Ei) = 0.106, for ao = 0.5, and Re{Ei) = 0.145, Re{E2) = 0.226, for ao = 2.25. 

The lifetime, tl of the quasibound states is determined from the imaginary part of the 
complex energy, Im{E), where the pole is found. Then 

rL = p (36) 

where F = 2Im{E) is the ionization rate. For the inverted Gaussian potential it has been 
found that with increasing strength of the driving field the ionization rate decreases in an 



oscillatory manner []TT],|12[. In Figs.(|3) and (|]) we show how the real and imaginary part of 
the quasibound state energies change as a function of ao, for ao ranging from to 6 a.u.. 
The incident particle can emit a photon and drop to a localized Floquet evanescent state. It 
is in this sense that in Fig. (^ we plot the real part of the quasibound state minus a photon 
energy and obtain results in agreement with those obtained in |lTT|,p!2|. 



B. Wigner-Smith delay times 

Wigner's |2D| one dimensional analysis on time delay in a quantum mechanical scatter- 



ing problem was generalized to multi-channel scattering by Smith |^ who introduced the 
Hermitian matrix 

Q^^S§ (37) 

and interpreted its diagonal elements Q„„ as the average delay experienced by a particle 
incident on the nth channel {S is the unitary scattering matrix). 

In what follows, we compute the diagonal elements of the Wigner-Smith delay matrix 
Q for the system currently under consideration. We first obtain the eigenvalues and eigen- 
vectors of the Floquet S'-matrix, see Eq . (|3T|) , which is a 2{np + 1) x 2{np + 1) matrix when 
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the system is truncated to + 1 propagating modes. The 2 (rip + 1) eigenvalues of the 
unitary Floquet S-matrix have unit magnitude and can thus be cast in the form e*^'^^\ 
where i = 1, 2(np + 1) and 6i{S) is the ith eigenphase as a function of the Floquet energy 
S. The eigenvector corresponding to the ith eigenvalue e^^' is denoted by \6i >. We note 
that the transmission coefficients r„ and Ttot,n (see section III. A) as well as the Wigner- 
Smith delay times r^^, defined in what follows, are a function of incident electron energy 
E G [0, {up + l)ti;). However, the eigenphases 9i and the eigenvectors \9i > of the Floquet 
S-matrix are a function of the Floquet energy S G [0, u) in the sense that if the Floquet 
energy is defined as the incident electron energy at a higher propagating channel one finds 
the exact same eigenvalues and eigenvectors. The Floquet S-matrix can be written as 

2(np+l) 

S= J2 \0^>e''' <0.\. (38) 

2=1 

Each of the eigenvectors \6i > can be expanded in terms of the propagating free electron 
waves \kn > that we have used to construct the Floquet S-matrix, where < x|/c„ >= e*'^"^. 
That is, 16*2 >= Y.'^o^ Pn,i\kn > where pn,i =< kn\0i > and Pn,i = |Pn,iP is the occupation 
probability of the \6i > eigenvector on the > propagating channel and n = 0, ...,np for 
the right propagating modes and n = Up + 1, ...,2np + 1 for the left propagating modes. 
The occupation probability of the \di > eigenvector on a mode incident from the right is the 
same with that of the corresponding mode incident from the left, that is, P„ j = Pn-\-np+i,i 
for n = 0,...,np. Finally, for each eigenvector 16*4 > the total occupation probability is 
normalized to 1, J2'^o~^ Pn,i = 1, where modes incident from the right and the left are taken 
into account. For modes only incident from the right the normalization for the occupation 
probability takes the form Y^n=oPn,i = 0.5. From Eqs. (^) and ([38| ) and the fact that the 
eigenvectors \6i > form a complete set (the Floquet S^-matrix is unitary) one can show that 

2(np+l) ,n 2(np+l) ,^ 

< A;„|Q|fc„>=C= Yl -:7^<kn\0^><0^\kn>^r:^s= E ^n.^^F" (39) 

i=l i=l '^^ 

The Wigner-Smith delay times r^^ are the average times an electron incident on the nth 
channel with energy E G [nu, {n + 1)uj) is delayed due to its interaction with the laterally 
oscillating time periodic potential in the K-H frame. The Wigner-Smith delay times for 
propagating modes incident from the right are the same with those incident from the left, 
that is, = r^+"''+^ forra = 0, Up. 



In Figs.(0), ( p!OD and ( [TI] ) we plot the Wigner-Smith delay times for equal to 0.5, 2.25 
and 5.25, respectively, for modes incident from the right. In Table I we compare the Wigner- 
Smith delay times r^^, obtained from the Floquet S-matrix with the lifetime tl obtained 
from the poles of the r„ transmission coefficient in the complex energy plane and find them 
to be of the same order of magnitude |^2[. At the transmission resonances, the incident 
electron wave gets trapped by the oscillating potential, populating the quasibound states of 
the system. The delay of the incident electron wave at the transmission resonances shows up 
as peaks when plotting the Wigner-Smith delay times as a function of the electron incident 
energy. Note that as Oq increases from 2.25 to 5.25 the Wigner-Smith delay time of the 
1st quasibound state increases, evidence of stabilization. Another interesting observation 
is that for small incident energy E the electron has positive Wigner-Smith delay times, for 
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ao = 0.5, but the electron has negative Wigner-Smith delay times for strong driving fields 
cto = 2.25 and 5.25. These last can arise physically either from reflection of the incident 
electron before it enters the scattering region or from its acceleration and swift passage 
through the negative potential [^. In addition, from Figs.(P), (^) and (0) we see, as 



expected, that when the electron is incident on higher Floquet channels it delays less and 
less until for very high energies it is not affected by the potential and the delay time is zero. 

Now, let us briefly comment on the truncation error of our numerical calculations. The 
number of Floquet channels N was chosen for each value of ao so that the error due to 
truncation remains small. The truncation error for the elements Rn' ^n/Tn' ,ni where n\n = 
0,...,np, is smaller for = and it increases as n approaches Up, where Up is the last 
propagating mode and thus the mode with the larger electron incident energy. Thus, the 
truncation error for the transmission coefficients T„ and Ttot,n computed in section III. A 
is smaller than the error for the Wigner-Smith delay times computed in this section. An 
estimate of the truncation error is given by 1 — iT^p ^^p as a function of incident electron 
energy. In all our calculations the truncation error is kept in the order of lO""^ for oq = 
0.5, 2.25 and 10"^ for ao = 5.25 so that our results are reliable. As we increase ao we need to 
consider a larger number of Floquet channels N to maintain a small truncation error in our 
numerical calculations making it computationally challenging to compute the Wigner-Smith 
delay times for large values of ao- 



C. Classical-Quantum Correspondence 

When we plot the eigenphases of the Floquet S-matrix as a function of the electron 
Floquet energy S we notice that the eigenphases undergo an increasing number of "avoided 
crossings" with increasing strength ao of the driving field, see Figs.(|T2|), (|14|), and (|1^). As we 



show in what follows, we believe that these "avoided crossings" are a quantum manifestation 
of chaos in the classical phase space. 

Let us explain what we mean by the term "avoided crossing" in terms of the occupation 
probabilities, Pn^i, (defined in section III.B) of the \6i > eigenvector on the > propagating 
channel. In what follows we consider the occupation probabilities only for modes incident 
from the right, that is, n = 0, Up. In Fig.(|l^) the eigenphases 6i and 62 undergo a repulsion 
when the Floquet energy is equal to the transmission resonance, S = 0.106, for ao = 0.5. 
For very small values of ao the eigenphases cross each other without repelling. It is only as 
we increase the strength of the driving field that the eigenphases undergo a repulsion which 
we refer to as an "avoided crossing". We describe quantitatively the "avoided crossing" 
between the eigenphases 61 and 62 in terms of occupation probabilities. In Fig. (13a) we plot 
the occupation probabilities Po,i) -^1,1 P2,i of the \9i > eigenvector on the propagating 
channels n=0,l,2 and in Fig. (13b) we plot the occupation probabilities Po,2, Pi, 2 and P2,2 
of the 1^2 > eigenvector on the propagating channels n=0,l,2 as a function of Floquet 
energy S. Before the "avoided crossing", S < 0.106, \6i > has support mainly on the second 
propagating channel, Pi^i ~ 0.5, and {62 > on the first propagating channel, Po,2 ~ 0.5, while 
after the "avoided crossing", S > 0.106, \9i > has support mainly on the first propagating 
channel, Po,i ~ 0.5, and {62 > on the second propagating channel. Pi, 2 ~ 0.5. This total 
exchange of character is what we refer to as a sharp "avoided crossing". Note that for 
ao = 0.5 the propagating channels involved in the "avoided crossing" are mainly n = 0, 1. 
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As the strength of the driving field ao is increased an increasing number of propagating 
channels undergo "avoided crossings" , as shown in Tables II and III. 

For increased strength of the driving field the number of avoided crossings increases, 
see Figs.(|l4D, (|l^) and it can be that more than two eigenphases participate in an "avoided 
crossing" for a certain Floquet energy. For example, this is the case for the "avoided crossing" 
at S ^ 0.14, for ao = 2.25, where there are three eigenphases 6i, 62 and 6q interfering, 
see Fig.([T3|). In Figs.(|T^ and ( PTBD we plot the eigenphases of the Floquet S-matrix as a 
function of S for ao = 1.25 and ao = 2.25, respectively, to show the increase in the number 
of "avoided crossings" with increasing strength of the driving field. In Tables II and III we 
present for ao = 1.25 and ao = 2.25, respectively, the eigephases which undergo "avoided 
crossings" for different Floquet energies S and the propagating channels n = 0,1,... for 
which the occupation probability Pn,i is substantial. In Tables II and III, the channels 
with small occupation probabilities are indicated as subscripts to the channels with large 
occupation probabilities. This is only an approximate picture but helps us visualize how 
the eigenphases change character at the "avoided crossings" . For example, from table III, 
we obtain an approximate picture how the eigenphases 61, 62 and 6q participate in the 
"avoided crossing" at £^ ~ 0.14. For S = 0.07, the eigenvector \6i > has support mainly on 
the propagating channel n = 1 and less on n = 0,2,3, the eigenvector {62 > has support 
mainly on the propagating channel n = and less on n = 1, the eigenvector \6q > has 
support mainly on the propagating channel n = 3 and less on = 5, 1. For S = 0.145, the 
eigenvector \6i > has support mainly on the propagating channel n = and less on n = 1, 
the eigenvector > has support mainly on the propagating channels n = 1, 3 and less on 
n = 5, the eigenvector \6q > has support mainly on the propagating channels n = 1, 3. Thus, 
there is an exchange of character among the eigenphases 61, 62 and 6q expressed in terms of 
the mainly interfering channels n = 0, 1, 3, 5 but it is not a complete exchange as in the case 
of the sharp "avoided crossing" at ao = 0.5, see Fig. (|l2D . The "avoided crossings" we have 
just described for the open quantum system under consideration, the inverted Gaussian in 
the presence of a driving field, are analogous to what was seen in a bounded chaotic system 
p3| where the authors also discuss two different types of "avoided crossings" . 

We now turn to the classical dynamics of the inverted Gaussian potential in the presence 
of the driving field. Figs.(0) are strobe plots of the phase space dynamics, for constant 
frequency u = 0.236 and increasing strength of the driving field, ao is equal to 0.5, 1.25 
and 2.25. The strobe plots are drawn by evolving a set of trajectories, with different initial 
conditions, and plotting the location of each trajectory at time intervals (t„ = 27rm/a;, 
m = 1, 2, ...) equal to the period of the driving field. We indicate the location of the period- 
1 periodic orbits with filled squares. The strobe plots are drawn in the Lab frame, see 
Eqs.([l|), (H) and (|^), and are exactly the same with those in the K-H frame except that in 



the Lab frame the x axis is shifted by ao p4|- driving field is present, the motion is 



regular and bounded for negative energies, while it is unbounded for positive energies. When 
the driving field is turned on, the KAM tori in the regular island around x = 0, p = start 
breaking up as ao is increased and chaotic motion sets in. For ao = 0.5, see Fig. (16a), the 
classical phase space is mixed. There are two islands around the two stable periodic orbits 
but there are also chaotic trajectories. As ao is further increased the remaining islands are 
very small, see Fig. (16b), until they totally disappear, see Fig. (16c), and the phase space 
in the scattering region becomes dominated by chaos. In addition, in Figs. (17a) and (17b) 
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where the initial values of the classical momenta are chosen to correspond to the middle 
of the Floquet propagating channels, we find that as is increased more trajectories get 
pulled into the chaotic region of the classical phase space. Correspondingly, in the quantum 
treatment of the scattering problem we have seen that as the strength of the driving field is 
increased the eigenphases of the Floquet S'-matrix undergo an increasing number of "avoided 
crossings" where more Floquet channels contribute to the scattering process. We thus believe 
that the "avoided crossings" are a quantum manifestation of the breaking of the constants 
of motion and chaos setting in in the classical phase space. 



IV. CONCLUSIONS 

In this paper, we have studied the scattering of electron waves from an inverted Gaussian 
potential, used to model the atomic potential, in the presence of strong time periodic electric 
fields. Using Floquet theory, we have constructed the Floquet S-matrix in the K-H frame, 
where asymptotic states can be defined. We have computed the transmission resonances, 
for different strengths of the driving field, and shown that they are associated with zero-pole 
pairs of the Floquet S-matrix in the complex energy plane. We have also computed the 
Wigner-Smith delay times which is a different way to "probe" the complex spectrum of the 
open quantum system. Finally, we have shown that the eigenphases of the open quantum 
system undergo a number of "avoided crossings" as a function of the electron Floquet energy, 
that increases with increasing strength of the driving field. We believe that the "avoided 
crossings" are quantum manifestations of the destruction of the KAM surfaces and the onset 
of chaos in the classical phase space. 

Acknowledgement: We wish to thank the Welch Foundation, Grant No F-1051, NSF 
Grant INT-9602971, and DOE contract No. DE-FG03-94ER14405 for partial support of this 
work. We also thank the University of Texas at Austin High Performance Computing Center 
for use of their computer facilities. Finally the authors thank A. Gursoy, T. Timberlake and 
A. Shaji for helpful discussions. 

Appendix A 



In what follows starting from Eqs.f^Tl), (|28|) , ( P9[) and (^) we obtain the outgoing proba- 
bility amplitudes of the propagating modes in terms of the incoming probability amplitudes 
of the propagating modes. 

Using = K and K '^K' = -K^ , we eliminate from Eqs.(g3) and (^) and 
obtain 

= M-^L^C + M^L^D, (40) 

where = K_X^_l^ - K'^'^X^^^'^ and L2 = K_X^l^ - K'^'^X^"^^'^. Using R-^ = K_ 
and K'j^ ^K'_ = —Kt, we eliminate from Eqs.(^) and (|30|) and obtain 

2KlJA''' = M-^L^C + M-^L^D, (41) 
where X3 = K_X^1^ + K'^'^X'^^^'^ and = K_X^^^ + K'^'^ X^'^'^'_. From Eqs.(|OD and 



(PI) we express C, D in terms of A*" and S*" as follows 
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C = 2G[L^^K'iJJ\fB''' - Ll^K'tJMA% (42) 

D = 2H[L^^KIjMB''' - L^^KlJAfA% (43) 

where G = [L^ ^Li - Ll^L^^-^ and H = [L^^L2 - L^^L^]'^. Substituting Eqs.(||) and 
(H) in Eqs.(|7D and (|2|) yields A°"* and B"""* in terms of A*" and 

^out ^ j^-^tWA'"" + Af-^rAfB'", (44) 



or equivalently 



where 



ATA""* = rWA''' + tAfB'"" 

ATB""* = t'MA"' + rATs'", (45) 



r 



+ 2K_X^l^GL2^Kl + 2K_X^^^HL^^Kl]J 



t' = -2[K^X^l^GL-^^Kl + K xf HL7^^Kl]J 

t = 2[K^X^}^GL^^Kl + K X^^ HL^^Kl]J (46) 

From Eqs.(^6D, due to the multiplication on the right by the N x N matrix J, we find that 
the N X N matrices r', r, t' and t are of the following form 

\ "pe ' PP / \ "pe ' PP / \ "pe ''pp / \ "pe •'pp / 

where the matrices r'^p, rep, t'^p and tep have dimensions rze xn^ + l and the matrices r'pp, rpp, 
t'pp and tpp have dimensions + 1 x + 1, respectively, (ng is the number of the evanescent 
modes and + 1 is the number of the propagating modes). The matrices Oge and Ope have 
dimensions rif. x and + 1 x ng, respectively, and they have zero elements because the 
amplitudes 6*" and a™ of the evanescent modes are zero, 6^" = a*„" = 0, for n = —rie, —1. 
In addition, the matrices Af, A*", A°"*, B*" and can be written as 

\r — i ee Ogp A \in _ { A^ \ \out _ { A^ \ j-,in _ ( \ j-fout _ ( -^e i 

(48) 

where the matrices Afee-, J^pp have dimensions rie x Ue and rip + 1 x rip + 1, respectively. 
The elements of the rie x 1 matrices A*", A""**, and B""* are the amplitudes of the 
evanescent modes. The elements of the np + 1 x 1 matrices A^", A°"*, and are the 
amplitudes of the propagating modes. Using Eqs. (|47|) and (|48|) we write Eqs.(|45|) as follows 



Kf Aout _ I Kj Ain,f Kf r>in 
•/V ee-'^e — ' ep"^* PP-^p ' ''ep-'V pp-'-'p 

Kf AOut _ t Kf Ain _i_f\f p>in 
j\ pp-m-p I ppJ\ pp-'^p ' '^pp^^ pp-'-'p 

MeeBT' = KpATppA;- + VepMppB^ 



MppB^' = t'ppAfppA;- + VppAfppBir (49) 
From Eqs. (^91) we obtain the Floquet S-matrix given in Eq . (^) . 
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Appendix B 



In sections II. B, C and D we have formally constructed the Floquet S-matrix in terms 



of the functions Xnjix), with n = — ng, 0, rZp and j = 1,...,2N, which are linearly 
independent functions in the scattering region II. For the inverted Gaussian potential the 
functions Xn,j{x) can only be obtained numerically. In Appendix A we formally expressed 
the matrices r and t in terms of the functions Xn.ji^)- Numerically, though, it is not efficient 

In what follows we outline the numerical method [01 we 



to compute the functions x 
use to obtain the N x N matrices r and t for electron waves incident from the right. 
The wavefunction in the asymptotic regions I and III is given by Eqs. ([T8|) and ( |T9| 



with 



a^j^ = 0, since we only consider electron waves incident from the right. We can then write 



0n(-a;o) = < 

with 6^" = for n = —rig, — 1. Next, we write Eqs.(^) in matrix form using Eqs. 
(p5|) and (^61) as follows 



(50) 



(51) 



cPi-xo) = K{-x^)J\fA 



out 



where 



/e 



— rig ^ 



Kix) 



-ikgx 



y Q g iknpX j 

Using Eqs.dl) with A'" = 0, we write Eqs.(|l|), (||) as follows 



(52) 



(53) 



(54) 



0(-xo) = K{-xo)tAfB' 



(55) 



Next, we use Eq.(^) to numerically propagate 0(x) from x = —Xq up to x = Xq according to 
the Numerov algorithm [^. From Eq. (|55|) we see that in practice we numerically integrate 
the K{x) matrix from x = — xq up to x = xq, since t, Af and -B™ are constant matrices. Let 
us indicate by K{xq) the numerically integrated matrix K{x) at x = xq. Then, matching 
the wavefunction and its first derivative at x = Xq and using Eqs. (03) we obtain 



k{xo)tAfB' 
k\xo)tArB' 



K+rMB''' + K JMB' 
K'^rAfB''' + K' JMB' 



(56) 



or equivalently 
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k{xo)t = K^r + K_J 

K{xo)t = K'^r + K'_J. ^ ' 

From Eqs.(|57D we find 

r = [K{xo)-'K+ - k\xo)-'K'^]~'[K{xo)-'K^ - K{xor'K'_]J 

t = [K.K{xo) - K'^'^K{xo)]-'[Kl - K'^^K' ]J, ^ ' 

where we have used the relation = K^. The matrices r and t given by Eqs. (pH]) are 
of the form shown in Eg .(^71) and thus we can extract the matrices Vpp and tpp. Following 
section II. D, we then obtain the matrices R and T. Finally, using the symmetry property 
of the Floquet S-matrix given in Eq. (|3^) we find the matrices R' and T' for electron waves 
incident from the left. 
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LIST OF TABLES 



1. Table L The Wigner-Smith delay times r^^ compared to the lifetime ti for the 1st and 
2nd quasibound states for equal to 0.5, 2.25 and 5.25. 

2. Table II: For ao = 1.25 we retain n = —6, 0, 6 channels and obtain 14 eigenphases 
from the 14 x 14 Floquet S-matrix. In Table II we only show the five eigenphases 
participating in the "avoided crossings" at different Floquet energies £, see Fig, (p^ . 
For each of the five participating eigenphases 9i,...,9^ we display the propagating 
channels n = 0,1,... with substantial occupation probability P„,j. The propagating 
channels involved in the "avoided crossings" are = 0, 1, 2. 

3. Table III: For ao = 2.25 we retain n = —12, 0, 12 channels and obtain 25 eigen- 
phases from the 25 x 25 Floquet S-matrix. In Table III we only show the seven 
eigenphases participating in the "avoided crossings" at different Floquet energies £, 



see Fig. ([T5|) . For each of the seven participating eigenphases 6i, ...,9j we display the 
propagating channels n = 0,1,... with substantial occupation probability P„ j. The 
propagating channels involved in the "avoided crossings" are ?t, = 0,1,...,5. 
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1. Figurel: The fourier components '^^'^ (a.u.) of the inverted Gaussian potential as 
a function of the one space dimension x (a.u.) in the K-H frame, for = 2.25 a.u.. 

2. Figure2: Not drawn to scale, are shown in the K-H frame the asymptotic regions I x G 
[xo,+oo) (a.u.) and III x e {—oo,—xq] (a.u.), where the potential is asymptotically 
zero, and the scattering region II where the inverted Gaussian potential oscillates 
laterally. In regions I and III, we also show the Floquet channels, denoted by dotted 
lines, and the incoming and outgoing electron waves, denoted by solid arrows. 

3. FigureS: The transmission coefficients T„ and Ttot,n, respectively, as a function of 
electron incident energy E, with E G [0, 2uj) a.u., for ao = 0.5 a.u.. There is only one 
Fano transmission resonance at E = 0.106 a.u., associated with the first quasibound 
state. 

4. Figure4: The transmission coefficients T„ and Ttot,n, respectively, as a function of 
electron incident energy E, with E G [0,2a;) a.u., for ao = 2.25 a.u.. There are 
two Fano transmission resonances ai E = 0.142 a.u. and E = 0.225 a.u., associated 
with the first and second quasibound states, respectively. The second order Fano 
transmission resonances for E > u are more prominent than those for ao = 0.5 a.u.. 

5. Figure5: The transmission coefficients Tn and Ttot,n as a function of electron incident 
energy E, with E G [0,20;) a.u., for ao = 5.25 a.u.. There are two Fano transmission 
resonances at E = 0.185 a.u. and E = 0.219 a.u., associated with the first and second 
quasibound states, respectively. The second order Fano transmission resonances for 
E > uo are more prominent than those for ao = 2.25 a.u.. 

6. Figure6: Contour plot of the transmission coefficient r„ in the complex energy plane 
for ao = 2.25 a.u.. The dark— >-light areas correspond to increasing values of r„. There 
are two zero-pole pairs each associated with the Fano resonances in Fig.(^. From the 
poles we determine the real part and the lifetime, tl, of the first and second quasibound 
states. 

7. Figure?: Real part of the first (squares) and second (dots) quasibound states minus 
a photon energy as a function of ao (a.u.). The real part of the quasibound states is 
found from the poles of the transmission coefficient r„ in the complex energy plane. 

8. Figure8: Ionization rate of the first (squares) and second (dots) quasibound states as 
a function of ao (a.u.). The imaginary part of the quasibound states is found from the 
poles of the transmission coefficient r„ in the complex energy plane. 

9. Figure9: The Wigner-Smith delay times, r^^, as a function of electron incident energy 
E G [0, 20;) for ao = 0.5 a.u.. There is one peak at E = 0.106 a.u. and smaller peaks 
at higher order resonances, associated with the Fano resonance in Fig.(|^). For small 
incident energy, E, the Wigner-Smith delay time is positive. 



21 



10. FigurelO: The Wigner-Smith delay times, r^^, as a function of electron incident energy 
E G [0, 3uj) for = 2.25 a.u.. There are two peaks at E = 0.142 a.u. and E = 0.225 
a.u. and smaller peaks at higher order resonances, associated with the two Fano 
resonances at Fig.@. For small incident energy, E, the Wigner-Smith delay time is 
negative. 

11. Figurell: The Wigner-Smith delay times, r^^, as a function of electron incident energy 
E G [0,40;) for = 5.25 a.u.. There are two peaks at E = 0.185 a.u. and E = 0.219 
a.u. and smaller peaks at higher order resonances, associated with the two Fano 
resonances at Fig.(|^). For small incident energy, E, the Wigner-Smith delay time is 
negative. 

12. Figurel2: The eigenphases 9i (rad) as a function of Floquet energy S (a.u.) for ao = 0.5 
a.u.. For oq -C 0.5 a.u. the eigenphases 9i and 02 intersect each other as a function of 
Floquet energy £ (a.u.). It is only as ao is increased that the eigenphases repel as a 
function of S (a.u.) and form an "avoided crossing", indicated by an arrow. 

13. FigurelS: For ao = 0.5 we retain n = —6, 0, ...6 channels and obtain 14 eigenphases 
from the 14 x 14 Floquet S-matrix. Only two eigenphases 6i and 62 participate in the 
"avoided crossing". The eigenphases 9i and 6*2 exchange character completely at the 
sharp "avoided crossing" shown in Fig. ([12|) . To show, quantitatively, how the character 
exchange takes place we plot in a) the occupation probabilities Pq^i, Pi^i and P2,i of 
the 1^1 > eigenvector on the propagating channels n = 0, 1, 2 and in b) the occupation 
probabilities Po^2, Pi,2 and P2,2 of the \92 > eigenvector on the propagating channels 
n = 0, 1,2 as a function of Floquet energy S (a.u.). Before the avoided crossing the 
eigenvector \9i > has support on channel n = 1 and the eigenvector {62 > mainly 
on channel n = 0, while after the avoided crossing the eigenvector \9i > has support 
mainly on channel n = and the eigenvector {62 > on channel n = 1, thus exchanging 
character completely. 

14. Figurel4: The eigenphases 6i (rad) as a function of Floquet energy S (a.u.) for ao = 
1.25 a.u.. The eigenphases 9i,...,9^ participate in the "avoided crossings" shown in 
Table II. 

15. Figurel5: The eigenphases 6i (rad) as a function of Floquet energy £ (a.u.) for ao — 
2.25 a.u.. The eigenphases 6*1,. ..,6*7 participate in the "avoided crossings" shown in 
Table III. 

16. Figurel6: Strobe plots of the classical dynamics, for the inverted Gaussian in the 
presence of the driving field, in the Lab frame for a) ao = 0.5, b) ao = 1.25 and c) 
ao = 2.25. The initial conditions used to generate the plots lie on the line p = as 
well as on the lines with — 1 < p < 1. The location of the period-1 orbits are indicated 
by filled squares. The period-1 orbits are located at a) —1.87, —0.36 and 2.46 b) 2.84 
and c) 1.88, 3.38 and 4.54. For very small values of the driving field ao (not shown) 
there is a large regular island around the region at a; = 0, p = 0. As ao is increased 
to 0.5 there are two regular islands reduced in size indicating the destruction of the 
KAM tori. As ao is further increased to 1.25 and 2.25 the regular islands disappear 
and the phase space is dominated by chaos. 
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Figure 17: Strobe plots of the classical dynamics, for the inverted Gaussian in the 
presence of the driving field, in the Lab frame for a) ao = 0.5 and b) = 2.25. 
The initial conditions of the classical momenta used to generate the plots are chosen 
to correspond to the middle of the Floquet propagating channels, that is the initial 
conditions lie on the lines p — ±^^2(0. 13 + nuj) with n = 0,l,...,8. As the strength of 
the driving field is increased more trajectories get pulled in the chaotic region in the 
classical phase space. 
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